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1 INTRODUCTION 



Recently Ifjurkert fc Tremaind (|2010h (hereafter BT) have 
presented evidence that the total population of globular 
clusters in a galaxy is almost exactly proportional to the 
mass of its central supermassive black hole. In rough terms 
such a relationship may not, on the surface, be terribly sur- 
prising because (as BT discuss) Ngc is nearly linearly re- 
lated to the luminosity or total mass of the host galaxy (e.g. 
Harris 1991) , and M, scales nearly l i nearly with the galax y 
bulge mass jWehner fc Harris! |2006| ; iFerrarese et alj 12006). 
At a more detailed level, however, the behavior of the GC 
population with galaxy s ize and type is m uch more complex 
jBrodie fc Straderll2006l ; IPeng et al.ll2008l ). and it is not ob- 
vious that the Ngc -M. correlation should be as tight and 
strongly defined as it is - as BT discuss, more so than is 
warranted by the measurement uncertainties. A potentially 
more surprising result they find from their sample is that 
the total mass contained in the globular cluster system is 
nearly equal to that of the central black hole. 

In this paper, we investigate the Ngc vs. M, correlation 
further and discuss additional implications arising from it. 



* E-mail: glharris@astro.uwaterloo.ca; 
ris@physics. monaster. ca 



ABSTRACT 

We explore the relation between the total globular population in a galaxy (Ngc) and 
the mass of its central black hole (M,). Using a sample of 33 galaxies, twice as large 
as the original sample discussed by Burkert & Tremaine (2010), we find that Ngc for 
elliptical and spiral galaxies increases in almost precisely direct proportion to M t . 
The SO-type galaxies by contrast do not follow a clear trend, showing large scatter 
in M, at a given Ngc- After accounting for observational measurement uncertainty, 
we find that the mean relation defined by the E and S galaxies must also have an 
intrinsic or "cosmic" scatter of ±0.2 in either log Ngc or log M». The residuals from 
this correlation show no trend with globular cluster specific frequency. We suggest that 
these two types of galaxy subsystems (central black hole and globular cluster system) 
may be closely correlated because they both originated at high redshift during the main 
epoch of hierarchical merging, and both require extremely high-density conditions for 
formation. Lastly, we note that roughly 10% of the galaxies in our sample (one E, one 
S, two SO) deviate strongly from the main trend, all in the sense that their M, is at 
least 10 x smaller than would be predicted by the mean relation. 

Key words: globular clusters, black holes, galaxies 



2 THE DATABASE 

The BT result was based on just 16 galaxies: 11 ellipticals 
(E), 3 spirals (S), and 2 lenticulars (SO), and only 13 were 
used in deriving the numerical relation between total cluster 
population and central black hole mass. Given that the num- 
ber of galaxies with either measured Ngc or M, is consider- 
ably larger than this, we chose to try and add to the sample. 
We carried out an independent search of the literature using 
galaxie s included i n the the M. databases of Gultekin et al.l 
2009d ), lGrafiaml (|2008l ). and lGebhardl |2004l ). whose glob- 
ular cluster systems (GCSs) had also been studied and for 
which it was possible to make reasonable estimates of Ngc- 
As a result we now have a database of 33 galaxies: 21 E, 
4 S, and 8 SO, doubling the BT sample. Many other galax- 
ies have published measurements of either one or the other 
(more than 70 with M. and well over 200 for Ngc), but 
these two types of observational studies have been carried 
out originally for very different purposes, so it is perhaps 
not suprising that the overlap between them is as yet quite 
limited. 

In Table 1 we list the results of our compilation. Succes- 
sive columns give the galaxy name, galaxy type, absolute V 
magnitude, total cluster population Ngc and data source, 
specific frequency (number of clusters per unit galaxy lu- 
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minosity), black hole mass M. and source. References for 
the Ngc and M. values are provided in the footnote. For 
a few of the galaxies in this list, the original GC study 
quoted only a value for the number of directly observed 
GCs and not a total cluster population. For these we re- 
constructed a total, correcting for spatial and magnitude 
incom pleteness with standard method s described, for exam- 
ple, in lHarris &: van den Berghll|l98ll) . In those cases we ref- 
erence both the original GC study and this paper. In some 
of these cases the resulting Ngc value we believe to be ac- 
curate to only a factor of two, but these are still quite useful 
for the purposes of this discussion; and it should be noted 
that very similar uncertainties also apply to many of the 
black-hole masses. For a thorough recent re view of the cur- 
rent M , database and analysis issues, see iGultekin et al.l 
|2009bl ). 

The key plot of log Ngc vs. log M. for all galaxies 
in Table 1 is shown in Figure 1. If the two quantities are 
exactly proportional then the mean relation in this log-log 
graph will have unit slope, and the data do follow a trend 
very close to that. Thus the first conclusion to be drawn 
from this larger dataset is that we verify the essential result 
of BT, now extending over a factor of 10 3 in both quantities. 
The main spine of the proportionality is defined primarily 
by the E galaxies which make up more than half the sample. 

Though a trend clearly exists, it is also true that some 
galaxies lie well off the range of scatter defined by the ma- 
jority. Particularly disparate are the points for NGC 5128, 
the Milky Way (MW), and the SOs NGC 4382 and NGC 
7457. All have black hole masses which are at least an order 
of magnitude too small for their cluster populations when 
compared to the main trend. For the MW, both M, and 
Ngc are very well known, and we cannot arbitrarily choose 
to treat the cluster population in any different way from the 
other galaxies in the list (for example by counting only the 
"bulge" GCs as BT suggest). In fact, given that clusters in 
the central bulge are generally more difficult to detect in 
external galaxies, a more reasonable 'correction' to the MW 
sample would be to eliminate some fraction of the bulge clus- 
ters. But since these make up ~ 30% of the total population 
at most, any such correction would still leave the MW as an 
outlier in Figure 1. 

For NGC 5128, the cluster population estimates have 
consistently been in the rang e of 1000-2000 and t he most 
recent determination of 1300 (jliarris. G.L.H.|[2010h is now 
more soundly based than the GC totals for many of the 
other systems in Table 1. To force this galaxy back to the 
mean line, either Ngc would need to be smaller by a fac- 
tor of at least 3 (an impossibility since more than 600 nor- 
mal GCs have been in dividually found and confirmed; see 
IWoodlev et all l|2010bh : or the BH mass would have to be 
far larger (also now unlikely, as mentioned above). The dust 
lane has made determination of M, rather challenging, but 
it does seem that the mass u sed here is one of high confi- 
dence l|Cappellari et all 120091 1. Lastly, in the cases of NGC 
4382 and NGC 7457 the BH masses are solid upper limits 
l|Gebharddl20ld ) and their clusters would essentially all have 
to disappear to make these galaxies fall anywhere near the 
main line. 

Two other objects worth special comment are the com- 
pact ellipticals M32 and NGC 4486B which have vanishingly 
small GC populations. Both are satellites of much larger 
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Figure 1. log Ngc 1s plotted against log M. for all 33 galaxies 
in our sample. E galaxies are filled red triangles, S galaxies are 
filled blue circles, and open green triangles are SO galaxies. For 
three galaxies only upper M m limits are measured and these are 
shown with the horizontal arrows; for one we have only an upper 
limit on Ncc> shown as a vertical arrow. The five galaxies which 
deviate most strongly from the main trend are labeled; see further 
discussion in the text. 



systems (M31 and M87 respectively) and it has often been 
suggested that the GCs they might initially have had were 
removed by tidal stripping (for the outlying ones) or dynam- 
ical friction (for the ones in the core). For M32, we decided 
to make a fresh estimate of the number of cluster s it might 
have b y searching in the M31 cluster catalog of iLee et al.l 
(2008). There are 9 clusters within a projected distance of 
3 kpc of M32 and within ±100 km/s of its radial velocity. 
From the same database, to within the count statistics we 
find the same number of number of M31 clusters per unit 
area, within that velocity range and at the mean radial dis- 
tance of M32; therefore we conclude that M32 has no sta- 
tistically significant GC 'signal'. We adopt a conservative 
upper limit N G c < 10. For NGC 4486B, the GC popula- 
tion as measured in the comprehensive deep Virgo Cluster 
Survey jPeng et al . 2008) is indistinguishable from zero, and 
the BH mass measurement is also only an upper limit, so it 
is not clear whether or not this galaxy falls off the mean re- 
lation. In any case, we do not use these two small ellipticals 
in the following analysis. 



3 ANALYSIS 

Given that the BT results were based on a sample dominated 
by E galaxies, which make up more than half of our list, 
we decided to look at the data sorted by subtype before 
trying to find numerical correlations. In Figure 2 we plot log 
Ngc vs. log M. as before, but with E, S, and SO galaxies 
in different panels. A brief glance is enough to see that the 
distributions show certain differences for the three galaxy 
types: 
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Table 1. SMBH Masses and Globular Cluster Numbers 
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• The E galaxies span a range of three orders of magni- 
tude in both M. and Ngc, larger than for the SO or S types, 
and they show a strong correlation between log M. and log 
Ngc over that range, even in the low M, - Ngc part of the 
graph where there are only a few galaxies. 

• The four S galaxies cover a factor of only ~10 in 
Ngc and ~100 in M.. But, with the notable exception of 
the MW, the spirals appear to follow a trend similar to the 
ellipticals. 

• The SO galaxies are quite different, occupying a similar 
Ngc range to the spirals, but a much larger one in central 
black hole mass. No clear trend is present. Arbitrarily in- 
creasing the central black hole mass for the two SOs with 
logA/. <7 would place them closer to the main group but as 
already discuss ed above, the qu oted values of M, are upper 
limits for both llGebhardtl boiO). That SOs have rather small 



cluster populations (that is, low specific frequencies) relative 
to ellipticals has long been known and can be illustrated well 
by examining the homogeneous ACS Virgo Cluster Survey 
lists of Peng et al. (2008). Their sample contains 31 galaxies 
identified as SO: more than half have Ngc J; 100 and all but 
two are likely to have < 400 clusters. 



Because the E galaxies as a subclass define by far the 
most well determined correlation, we use them to define the 
baseline relation. We solve d for a best-fit relation using the 
X 2 estimator described by ITremaine et ail (2002) which ac- 
counts for the observational measurement uncertainties in 
both quantities for each galaxy individually. The 18 ellip- 
ticals (excluding only M32, NGC 4486B, and NGC 5128) 
yield a solution 
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Figure 2. log Ngc versus log Mm for the three types of galaxies (E, SO, S). 



logiV G c = (-5.78 ±0.85) + (1.02 ± 0.10)log— p- (1) 

Mq 

or, equivalently, 

A/rw 

log— p- = (-5.66 ±0.29) ± (0.98±0.10)logiV G c (2) 

Mq 

Including NGC 5128 does not change the best-fit relation 
significantly though it does increase the ^(min). 

In short, from this larger database we find (as did BT) 
that the mass of the central black hole scales nearly linearly 
with the globular cluster population in the galaxy. If we force 
the slope in Equation (1) to be exactly 1.00, then the best-fit 
relation would become 

Mbh/Mq = 4.07 x 10 5 N GC (3) 

or, expressed as a rule of thumb, about 250 globular clusters 
per 10 8 M Q black hole. 

Figure 3 shows both the best-fit relation of Equation (1) 
and the unit-slope relation, superimposed on the separate 
plots for the three types of host galaxies. Plotting both of 
these Uustrates that we have no compelling reasons to reject 
an exact 1:1 relation. It should be possible to improve the fit 
with the addition of more galaxies, particularly those with 
log M m < 8. 

The residual scatter about the mean (as indicated by 
the raw -^{min)) is noticeably larger than the random mea- 
surement errors would nominally permit. We have no partic- 
ular reason to assume that the various carefully done stud- 
ies in the literature have systematically underestimated the 
measurement uncertainties in either quantity, so we suggest 
that a real cosmic dispersion is also present. To bring the 
formal \ 2 P er degree of freedom down to 1.0 at the best-fit 
point, we find that we need to introduce an additional in- 
trinsic scatter of ±0.2 in either log M, or log Ngc, or in 
linear terms a factor of ~ 1.6. 

For a normal Gaussian-type globular cluster luminos- 
ity function such as in the Milky Way, with a mean 



at My = —7.3 and dispersion a(Mv) = 1-3 (e.g. 
Brodie & Strader 2006), the average cluster luminosity is 
L v = 1.2 x 10 5 L(T). With a mea n (M/L v ) = 2 (e.g. 
iMcLaughlin fc van der Marell (|2005l )) the mean GC mass is 
then 2.4 x 1O 5 M0, and the E-galaxy correlation becomes 
very roughly Mbh ~ 1.5MGc(tot). To well within a factor 
of two the total masses in these substructures are the same. 

The spiral galaxies, with the exception of the MW as 
discussed above, fall extremely close to the E galaxy line. 
In fact, adding the three 'good' spirals to the E-galaxy list 
and recalculating the solution leads to no difference from the 
relation given in Equation (1). The SOs are a different matter 
however; while a few fall nicely on or near the relationship, 
there are an equal number for which the M, - Ngc relation 
appears to be irrelevant. 

It has long b e en k nown (beginning with 
lHarris fc van den Berghl l|l98ll )) that the GC specific 
frequency Sn, i.e. the number of clusters per unit galaxy 
luminosity, varies with both g alaxy t ype and env ir onme nt 
(see also iBrodie fc Straded (|2006h ; iPeng et all (|2008h ). 
Since the galaxies in the list of Table I have a large range in 
Sn, we thought it would be worthwhile to ask whether any 
extra residual correlation of M, with Sn might be present. 
For example, consider two elliptical galaxies with the same 
total Ngc but different luminosities. By definition, the 
higher luminosity one must then have a lower Sn- But, we 
might also expect the higher-luminosity galaxy to have a 
higher mass M, if the central black hole mass is driven 
primarily by the host galaxy's mass. If this argument 
were valid, we should then see the opposite effect - i.e. 
an anticorrelation of M, with Sn- In Figure [3] we show 
the residuals in M, from the mean relation of Equation 
(1) plotted as a function of Sn- Interestingly no trend is 
apparent, or at the very least, any such trend is obscured 
by the measurement scatter in M, and the small number 
of points at high Sn- Of the two largest Sn galaxies, one 
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Figure 3. The same plots as in the previous figure, now with the best-fit relation for the E galaxies plotted in all three panels (red line). 
The black line in each panel is the relation with the forced assumption that the slope is unity. 



(M87) has a nearly "normal" M. and the other (NGC 1399) 
has a slightly lower than "normal" M, . In short, from the 
present data we see no extra correlation of central black 
hole mass with specific frequency. 

The most sparsely populated area of the correlation 
is at low mass. For M, much less than 1O 8 M0, the pres- 
ence of nuclear star clusters (NCs) may also become rele- 
vant |Graham fc Spitlerl [20091 ). Both BH and NC are part 
of t he more general class of Central Massive Objec ts (CMO; 
see IWehner fe Harris! 120061 ; iFerrarese et al.l 120061 ) and, as 
iGraham fc Spitlerl l|2009l )show. for some galaxies both are 
present at once. For intermediate-mass galaxies there is a 
gradual mass transition zone such that for M. <2x 1O 7 M0 
the NC mass tends to become dominant. Inspecting Table 
1 of Graham & Spitler, we note that for the MW and NGC 
7457 particularly, adding the NC mass to M. would bring 
them much closer to the mean E/S relation. For several other 
galaxies, however (NGC 1023, 1399, 3115, 4697), no impor- 
tant change would result. Additional data may help clarify 
whether or not the NC mass is an essential part of the pic- 
ture. 

As said earlier, the basic fact that the GC populations 
and central black hole mass should be correlated is not sur- 
prising, because in rough terms bigger galaxies should have 
bigger subsystems of all types. What is more intriguing is 
that for the E and S galaxies especially, the correlation 
is so tight and so nearly linear. BT suggest that mergers 
may be responsible for the growth of both these subcom- 
ponents of the galaxy. However, in a more general sense 
mergers will promote the growth of all subcomponents of 
the galaxy. We speculate that a more relevant link between 
black holes and globular clusters may be their age. Both 
the central BH and the majority of the GCs had their ori- 
gins at high and nearly similar redshift, during the major 
stages of hierarchical merging. Most GCs have ages in the 



range 10 — 13 Gyr, corresponding to redshifts z ~ 2 — 7 
jMarin-Franch et al.ll2009l;|Puzia et al.ll2005l ; IWoodlev et al.l 



l2010al ; iMuratov fc Gnedinl |2010T| . For large galaxies the 
seeds of the central black holes are in place before z ~ 6 
and grow by further gas accretion til l z ~ 2 and later (e.g. 
Alexander et al. 2005; Haiman et al. 2007; Volontcri et al. 
20031; iTanaka fc Haimanl 120091; IVolonteri fc Rees! 120051; 
Willott et alJ l20ld; iHeckman et al.1 12004|; ISilverman et al l 
hoOSl ; iMerloni et al.l l2010l : iKellv et alj |2010L among many 
others). Both types of structures clearly require formation 
conditions of extremely high gas density (> lO 5 M0/pc 3 
within small, parsec-sized volumes), which would have been 
easiest to generate in large numbers at very early times when 
gas cloud collisions and hierarchical mergin g events were fre- 
quent and energet ic (cf. the discussions bv lPeng et al.ll2008l ; 
iRhode et ail 120051 ). In this scenario, Nqc and M, should 
be closely correlated simply because they are both byprod- 
ucts of similar extreme conditions at high-density locations 
during the main period of galaxy formation. 

A less certain issue would be whether or not there is 
a direct causal link between them. For example, the early 
energy output from the dense, violent star formation condi- 
tions in the central regions of a large galaxy that might 
have accompanied the buildup of the massive BH might 
also have helped stimulate the formation of dense, mas- 



sive star clusters throughout the bulge and halo regions 
(e.g. lHarris et al.|[l998; Taniguchi et al.lll999l; rKrausell2002]; 



ISpringel et"al1l2005l ; iKauffmann et al.ll2007l ) . 

The rough equality of total mass in the two subsystems 
as shown above may, however, be largely a coincidence. The 
reason is that it applies only to the present stage of evo- 
lution of both subsystems. The total mass in the globular 
cluster system has been continually decreasing since their 
formation epoch: during the first several ~ 10 7 years of a 
young massive star cluster's life, it sheds a high fraction of 
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its initial protocluster gas due to SNe and stellar winds, and 
the highest-mass stars evolve and disappear. Over the subse- 
quent Hubble time, tidal stripping and evaporation remove 
most low-mass cluste rs, and the high er-mass ones are sig- 
nificantly eroded (see IVesperini|[2010l . for a review). These 
factors combined mean that the initial mass in the GC sys- 
tem should have been many times larger than it is now. By 
contrast, the BH mass has been continually growing over the 
same time through gas accretion and stellar capture, though 
its most rapid growth stage would have been at early times 
(cf. the references cited above). In other words, in their early 
stages the total mass in the GC population must have been 
much larger than the initial BH mass (an order of magni- 
tude or more). It is only now, after a Hubble time's worth 
of evolution in opposite senses, that their total masses are 
nearly crossing. 

The case of the SO galaxies is extremely puzzling. There 
seems to be no obvious way in which the data for these sys- 
tems is more uncertain than for the E and S galaxies, leaving 
us in the position of questioning the sample and/or classi- 
fication. As already discussed, the two most discordant SOs 
have black hole masses which are upper limits and there 
appears to be no obvious way in which they can be moved 
to more closely agree with the majority of the galaxies in 
our sample. Interestingly, they deviate from the norm in the 
same sense as the MW and NGC 5128 (discussed above). 
The four galaxies which are the most discordant then (1 
E, 1, S and 2 SO - roughly 10% of our sample) all devi- 
ate in the sense that their M, values are too small com- 
pared with Nqc', and without the two discordant SOs, the 
SO sample overall becomes less peculiar. However, all four 
"discordant" galaxies have normal specific frequencies for 
their galaxy type - i.e. Ngc is normal when compared to 
the luminosity /mass of the whole galaxy. This suggests to 
us that the formation of GCs is more strongly correlated 
with overall galaxy properties, and that in some cases M. is 
a later or less rapid development. 
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Figure 4. Residuals in M. from the mean relation in Equation 
(1) are plotted versus globular cluster specific frequency. E galax- 
ies are shown as solid red dots, S galaxies as blue triangles, and 
SO galaxies as green triangles. 

On observational grounds there is much room for im- 
provement to explore this intriguing correlation further with 
bigger samples. Galaxies with already-known BH masses can 
be imaged with deep wide-field photometry to determine 
their cluster populations, while the cores of many nearby 
galaxies with well determined Nqc values can be investi- 
gated to determine new BH masses. Some dozens of new 
datapoints can be added to the current list. Larger sam- 
ples should also allow exploration of wider ranges of galaxy 
environment and host galaxy type, and perhaps allow us to 
uncover more direct causal links between these old substruc- 
tures. 



4 CONCLUSIONS 

(i) We have investigated the correlation between central 
black hole mass in a gala xy and the total globular clus- 
ter population suggested bv lBurkert &: Tremaind (|2010h . By 
doubling the sample of galaxies and reconstructing the num- 
bers for Ngc and M, from the primary literature, we find 
results that confirm the basic trend they presented. 

(ii) Over three orders of magnitude, the elliptical and spi- 
ral galaxies in the sample define a relation whereby Ngc ~ 

i a nearly linear proportionality. By contrast, the 
SO galaxies in the sample do not follow the same trend, show- 
ing much larger scatter in M.. 

(iii) Of 33 galaxies in the current sample, three or four 
(i.e. ~ 10%) fall far enough from the mean relation to be 
considered genuinely anomalous. The Milky Way is one of 
these, with a central BH mass far too small for its cluster 
population. 

(iv) We suggest that the source of this correlation is con- 
nected with the epoch of origin of these structures. Both the 
central BH and the bulk of the GC population had their ori- 
gins at high redshift (z ~ 2 — 7), both requiring extremely 
high gas density conditions of lO 5 M0pc -3 and above. 
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